Bloch oscillations in lattice potentials with controlled aperiodicity 
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We numerically investigate the damping of Bloch oscillations in a one-dimensional lattice poten- 
tial whose translational symmetry is broken in a systematic manner, either by making the potential 
bichromatic or by introducing scatterers at distinct lattice sites. We find that the damping strongly 
depends on the ratio of lattice constants in the bichromatic potential and that even a small concen- 
tration of scatterers can lead to strong damping. Moreover, collisional interparticle interactions are 
able to counteract aperiodicity-induced damping of Bloch oscillations. The discussed effects should 
readily be observable for ultracold atoms in optical lattices. 
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I. INTRODUCTION 

The oscillatory motion of particles in a periodic po- 
tential, when subject to an external force, was predicted 
by Felix Bloch in 1928 1 . Bloch oscillations were first 
observed in the 1990s both in semiconductor superlat- 
tices [5] and in systems of laser-cooled atoms in optical 
lattices [31 0] . Since then, they have also been studied 
in atomic quantum gases [SHE] and have found appli- 
cations in ultracold atomic-physics-based precision mea- 
surements [9"HT2"] . 

Bloch oscillations are absent in solids due to fast damp- 
ing from scattering by defects and phonons. Their obser- 
vation in long-period semiconductor superlattices relies 
on oscillation periods that are shorter than the charac- 
teristic scattering lifetime, and even in those systems, 
damping of Bloch oscillations due to disorder is ubiq- 
uitous |131 114 ) . In contrast, such damping is absent in 
optical lattice systems which are inherently defect-free, 
allowing for the observation of a large number of oscil- 
lations [BHH]- Damping can however be induced by the 
introduction of disorder into the lattice potential, which 
can in principle be done with various techniques [TSHT§] . 
Also, in the case of quantum g damping of Bloch 

oscillations arises from mean-field interactions between 
weakly interacting, Bose-condensed atoms [7J [51 l2"Drl2"4"] . 
In the presence of both interactions and disorder, a re- 
duction of disorder-induced damping due to screening of 
disorder by the mean field has been predicted [25] . Ex- 
periments with ultracold atoms thus not only constitute 
a versatile testbed for behavior expected in solid-state 
systems but may also display novel effects. 

Recently, the damping of Bloch oscillations in a Bosc- 
Einstein condensate j^Hj has been observed for an optical 
lattice with a superimposed randomly corrugated opti- 
cal field. A related theoretical investigation of disorder- 
induced damping of Bloch oscillations has considered 
the case of Gaussian spatial noise [25]. In this paper, 
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we numerically investigate the dynamics of an atomic 
wave packet in a potential with a systematically degraded 
translational symmetry, considering two scenarios: The 
first is based on the use of a weak bichromatic poten- 
tial [H] of a variable wavelength ratio, and the second 
considers scatterers (impurities) pinned at single sites of 
the potential [17] . We find that the damping strongly 
depends on the ratio of lattice constants in the bichro- 
matic case and that even a small concentration of scat- 
terers can lead to strong damping. We also include ef- 
fects of the mean-field interaction and find that the rate 
at which damping of the Bloch oscillations occurs is re- 
duced, similar to the case of Gaussian disorder [25J. Both 
effects should be observable experimentally with existing 
ultracold-atom technology. 

This paper is organized as follows: After a brief discus- 
sion of fundamental aspects of tilted lattices in Sec. [TTJ 
Section |III| investigates the influence of aperiodicity on 
the damping dynamics of Bloch oscillations. Sec. [IV] ad- 
dresses the interplay of aperiodicity and the interaction 
between atoms. Conclusions are given in Sec. [Vj 

II. TILTED PERIODIC POTENTIALS 

The Hamiltonian for the motion of a particle in a one- 
dimensional periodic potential V(x) — V(x + a) with lat- 
tice constant a possesses a complete set of eigenfunctions 
that obey Bloch's theorem (p(x + a) = e ika ip(x). The 
corresponding eigenvalues E n {k) are periodic in momen- 
tum space and form energy bands, E n (k + K) — E n (k) 
(where n is the band index, hk the quasimomentum, 
and K = 2ir/a the width of the first Brillouin zone). 
Under the influence of an externally applied constant 
force ("tilt") F, the quasimomentum evolves as hk(t) = 
hko + Ft. Due to the periodicity of the energy bands, 
this results in oscillations of the particle's group veloc- 
ity Vg.n(k) — (1/h) dE n (k)/dk. These oscillations occur 
with a period Tg = 2tt/ujb = h/(aF) and a maximum 
displacement 2Ab,u — ^-n/F in coordinate space, where 
A„ = \E n (K/2) - E n (0)\ is the width of the n-th band. 

In the following, we consider particles that are con- 
fined to the lowest Bloch band, n=l, corresponding to 



FIG. 1. Evolution of the momentum-space density \i/j(k,t)\ 2 
(a) in a tilted periodic potential (b) in the presence of an 
additional periodic potential, and (c) in the presence of two 
localized scatterers as discussed in the text. Additional mo- 
mentum components emerge in (b) and (c), broadening the 
momentum distribution. The parameters in (b) are 7 — 0.01 
and a = l/v5; the scatterers in (c) are spaced 11 sites apart. 
A detailed explanation is given in the text. The density of 
\ip(k,t)\ 2 is normalized to 1; a corresponding color scale is 
shown on the right. 
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FIG. 2. Damped Bloch oscillations in a bichromatic lattice 
(7 = 0.005, a = 0.4). The collapse of the oscillations in 
coordinate space (a) is accompanied by breathing-mode exci- 
tations of the wave packet (b), and vice versa (see also text). 
The dashed line in (a) is the envelope Aexp(— rjt 2 ) + B. 



A. Bichromatic potentials 



sufficiently deep potentials and small enough tilts such 
that Zener tunneling [57] at the band edges is negligible. 
Using the split-operator method 28 , we perform numer- 
ical simulations of the dynamics of Bloch oscillations of 
a Gaussian wave packet 



^ = ° )= (2^)1/4 ex P 



(1) 



that evolves according to the Hamiltonian 



2m 



H 



+ V Q cos(Kx) + V(x) + Fx (2) 



for motion in a tilted periodic potential V(x) — 
Vq cos(Kx) + Fx that is modified by a weak additional 
potential V(x). 



III. EFFECTS OF MODIFIED PERIODICITY 



This section investigates the influence of the additional 
potential V(x), which either is a weak additional peri- 
odic potential with variable lattice constant, or arises 
from the local interaction with scatterers pinned at sin- 
gle sites of the tilted periodic potential. In the absence 
of V(x), the energies of neighboring sites differ by a fixed 
amount Fa — Hlub , corresponding to a spatially homoge- 
neous phase difference A0(f) = uist. In the presence of 
V(x), A(j) generally becomes position-dependent, leading 
to global dephasing and thus to a broadening of the wave 
packet in momentum space, as illustrated in Fig. [T] 



A tunable bichromatic potential is generated by the 
addition of 



V(x) = jVq cos(aKx), 



(3) 



with variable relative amplitude 7 <C 1 and lattice- 
constant ratio a. If a is a rational number a = p/q, 
the total potential has a periodicity A — aq. If further- 
more A exceeds the spatial range probed by the wave 
packet, the potential can be considered disordered. 

The evolution of the wave packet (Eq. [IJ typically ex- 
hibits collapses and revivals of center-of-mass oscillations 
that are coupled to breathing-mode excitations at twice 
the Bloch frequency as is shown in Fig. [2] (visible after 
a transient phase immediately following the switch-on of 
V) . The presence of a collapse and a revival of the Bloch 
oscillations is a consequence of the absence of dissipation 
in our model. To characterize the decay of Bloch oscil- 
lations, the numerical data for their initial decay, for a 
given 7 and a, are fitted with the function 



/(t) = Acxp(-^ 2 ) cos (u) B t) + B . 



(4) 



Results for the Gaussian-decay constant 77 are shown in 
Fig.0 

Clearly an increase in the perturbation amplitude 7 
leads to an overall increase in the damping of the oscil- 
lations. However, the dependence of the damping on the 
ratio a is less trivial. For integer values of a, the poten- 
tial retains its original periodicity and no dephasing of 
Bloch oscillations occurs. It is also suppressed for half- 
integer values of a (i.e. for q = 2), which correspond to 
a doubling of the lattice constant, and correspondingly a 
halving of the Brillouin zone, and of the corresponding 
Bloch period. For this case, the dynamics of the wave 
packet in the modified band structure can easily be visu- 
alized. The shape of the band is essentially that of V(x), 
but it is folded back into the new Brillouin zone, thus 
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FIG. 3. Damping constant 77 as a function of the lattice- 
constant ratio a, for three different depth ratios: 7 = 0.005 
(squares), 7 = 0.01 (diamonds) and 7 = 0.013 (circles). The 
inset shows the behavior of the decay constant in the vicinity 
of a = 1.5. 
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FIG. 4. Wave packet in a bichromatic potential with a = 
1/2,7 = 0.005. (a) Overlap | {ip(x,t = 0)|V>(M)) I, exhibiting 
the original periodicity of Tb, with small, growing, contri- 
butions at odd multiples of Tg /2. (b) Band structure of the 
potential. The inset shows the small gap at the Brillouin zone 
boundaries, with a calculated width A = 9.5 x 10~ 5 E r . The 
arrows with the corresponding numbers indicate the motion 
of the particle in k space. The wave function completes one 
Bloch cycle after Tb, going from 1 — s> 2 — > 3 — > 4, with most 
of the wave function tunneling across the gap. 



forming a closed loop, with a small splitting A a 7 at 
the new zone boundaries (cf. Fig. In one Bloch cy- 
cle, most of the wave function tunnels through the tiny 
gap from the lower to the upper portion at one boundary 
and then back to the lower portion at the other bound- 
ary (Bloch-Zener tunneling [3U]). The time needed 
for one such cycle is Tb, that is the Bloch period for the 
unperturbed potential V(x). This behavior can directly 
be seen in the time dependence of the overlap of the wave 
function with the initial wave packet | (ip(x,0)\ip(x,t)) |, 
which has a periodicity of Tg ', contributions at odd mul- 
tiples of Tg/2 only grow very slowly over a large number 
of cycles. 



The dynamics of Bloch oscillations in Figs. [2] and [3] are 
for a wave packet with width 00 = 20a in a lattice with 
depth V = 10 E r [where E r = (hK/2) 2 /2m is the recoil 
energy], and a small external force F — 0.011 Vol a. In 
coordinate space the amplitude of the Bloch oscillations 
is Ab = 0.34 a. The additional lattice V(x) is chosen 
to have a small relative amplitude 7 = 0.005, which is 
sufficient to cause noticeable damping already after a few 
Bloch cycles, depending on the lattice constant ratio a. 

Experimentally, such a bichromatic potential is 
straightforward to realize in the context of optical lat- 
tices, using two laser beams of different wavelength. With 
present-day tunable-laser technology, a large fraction of 
the range of a shown in Fig. [3] can be accessed. 



B. Scatterers on distinct sites 

For scatterers pinned at a set of distinct sites {n} of the 
lattice, we model the potential V as a sum of Gaussians, 



I/(x) = ^iexp[- 
{»} 



{x - x n ) 2 /{2d 2 )} 



(5) 



For the simulation, the parameters of the optical lattice 
and the tilt are chosen such that the size a of the wave 
packet and its oscillation amplitude Ab in the unper- 
turbed titled potential each cover a large number of lat- 
tice sites. A small number of scatterers are then ran- 
domly placed within the range of the wave packet's mo- 
tion. The amplitude A and width a of the scatterers are 
chosen such that the valleys of the unperturbed poten- 
tial Vo cos(Kx) are effectively filled up where the scat- 
terers are located. Their effect on the Bloch oscillations 
is shown in Fig. [5] The general trend with an increas- 
ing number of scatterers is a more rapid damping of the 
Bloch oscillations, with the details of the dynamics de- 
pending on their spatial arrangement. Clearly, already a 
small number of scatterers can lead to a rapid damping 
of Bloch oscillations. 

The results shown in Fig. [5] are for a comparatively 
shallow lattice (Vq = 1.4 E r ) and a tilting force F = 
0.022 E r /a, resulting in a large amplitude Ab = 16 a 
for undamped Bloch oscillations. This condition for the 
scatterers is well fulfilled by setting a = 7ra/20 and A — 
Er. 

In the context of ultracold atoms, the placement of 
scatterers pinned at single lattice sites can be achieved, 
for example by using atoms with two internal states in 
conjunction with a state-dependent lattice depth |31) . or 
by using two atomic species in a species-dependent opti- 
cal lattice 132 . 



IV. EFFECT OF INTERACTIONS 

The discussion so far has been restricted to noninter- 
acting particles. We now consider the case of a Bose- 
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FIG. 5. (a) Evolution of the packet position (x(t)) in the 
presence of randomly distributed scatterers. Increasing the 
number of scatterers (as indicated) increases the damping of 
Bloch oscillations. Each curve in a plot represents a differ- 
ent spatial configurations, (b) Overlap with the initial wave 
function for the case of randomly distributed scatterers. The 
parameters for (a) and (b) are given in the text. 



damping constant rj 7 up to a characteristic value g c , be- 
yond which an increase in g leads to an increase in rj. 
We interpret the reduction of 77 as a partial screening 
of the potential corrugations by the mean field [55] as g 
increases, which is eventually overcompensated by mean- 
field-induced dephasing [21"]l22]. 
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FIG. 6. Interplay between aperiodicity and mean-field inter- 
action, showing (a) evolution for g = (solid), g = 0.1 E r a 
(dotted) and g = 0.4 E r a (dashed), and (b) dependence 
of g on the coupling constant g, reaching a minimum at 
g c = 0.33 E r a 



The parameters for the simulation are the same as 
those in Sec. 3.1, leading to an associated coupling 
strength g c = 0.33 E r a. This value should be compared 
with an effective one-dimensional interaction parameter 
9 = 53£>/(27ra 2 L ) for a trapped atomic Bose-Einstein con- 
densate, where g^n = 4:irh 2 a s N/m (with atomic s-wave 
scattering length a s and atomic mass m), and a± ~ R, 
where R is the Thomas-Fermi radius. Already for a small 
condensate of 87 Rb atoms with N = 1 x 10 4 atoms and 
R = 5.3 /mi in an isotropic 50 Hz trap [33], and an op- 
tical lattice with a = 532 nm, we obtain g ~ 0.42 E r a, 
which is in the vicinity of g c . Hence, a significant mod- 
ification of the damping rate in a bichromatic potential 
due to mean-field effects can be expected; the coupling 
g depends, for example on the atom number N in the 
condensate, which is variable. Alternatively, an inves- 
tigation of the interplay is possible for species in which 
the mean-field interaction can be tuned via a Feshbach 
resonance. In this context we mention that Bloch oscil- 
lations with widely controllable mean-field interactions 
in a (monochromatic) optical lattice have recently been 
demonstrated with cesium condensates [71 |S] • 



Einstein condensate of N atoms with repulsive interpar- 
ticle interaction U(xi —Xj) = gS(xi — Xj). The evolution 
of the condensate wave function ip(x,t) = \/N 4>(x,t) 
in the bichromatic potential is then determined by the 
mean-field Hamiltonian 

H = Hq + -yV cos(aKx) + gN\<j>{x, t)\ 2 . (6) 

The simulation results (Fig. [6]) show that an increase 
in the coupling strength g leads to a reduction of the 



V. CONCLUSIONS 

We have numerically investigated the damping of 
Bloch oscillations resulting from a controlled breakdown 
of the periodicity of the lattice potential. The effects dis- 
cussed here, including the effects of the mean-field inter- 
action on disorder-induced dephasing, should be readily 
observable in experiments with ultracold atoms in optical 
lattices. 
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